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I. INTRODUCTION 

Recent advances in gravitational and quantum physics 
indicate that in order to reconcile the two fields with 
each other, a dramatic conceptual shift is required in 
our understanding of spacetime. In particular, the no- 
tion of spacetime as a continuum may need revision at 
scales where gravitational and electro-weak interactions 
become comparable in strength [l[ . For this reason there 
has been a recent revival of interest in approximating 
the spacetime with discrete coarse-grained structures at 
small, typically Planckian, length scales. Such structures 
are inherent in many models of quantum-gravity jSuch as 
spacetime foam , loop quantum gravity [1, [j, Q , non- 
commutative geometry 0, [7,, 8, 9i\ , black- hole physics [l0| 
or cosmic cellular automata lllTl^. fl^ fisj . 

Despite a vast gap between the Planck length (£p « 
1.6 • 10~^^ m) and smallest length scales that can be 
probed with particle accelerators (« 10"^^ m), the is- 
sue of Planckian physics might not be so speculative as 
it seems. In fact, probes such as Planck Surveyor [l^ 
or the related IceCube [13] — which just started or are 
planned to start in the near future, are supposed to set 
various important limits on prospective models of the 
Planckian world. 

One of the simplest toy-model systems for Planckian 
physics is undoubtedly a discrete lattice. Discrete lat- 
tices are routinely used, for instance, in computational 
quantum field theory [isl. [l9l . [20| . but with a few no- 
table exceptions [2l|, |22|. |23||. they mainly serve as nu- 
merical regulators of ultraviolet divergences. Indeed, a 
major point of renormalized theories is precisely to ex- 
tract lattice-independent data from numerical computa- 



tions. One may, however, investigate the consequences 
of taking the lattice no longer as a mere computational 
device, but as a bona-fide discrete network, whose links 
define the only possible propagation directions for sig- 
nals carrying the interactions between fields sitting on 
the nodes of the network. 

Recently one of us proposed a model of a discrete, 
crystal-like universe — "world crystal" [2^ [23 |. 
There, the geometry of Einstein and Einstein-Cartan 
spaces can be considered as being a manifestation of the 
defect structure of a crystal whose lattice spacing is of 
the order of £p. Curvature is due to rotational defects, 
torsion due to translational defects. The elastic deforma- 
tions do not alter the defect structure, i.e., the geometry 
is invariant under elastic deformations. If one assumes 
these to be controlled by a second-gradient elastic action, 
the forces between local rotational defects, i.e., between 
curvature singularities, are the same as in Einstein's the- 



ory [25|. Moreover, the elastic fluctuations of the dis- 



* Electronic address: 'jizba@physik.fu -berlin.de| 
t Electronic address: kleinert@physik. fu-berlin.de| 
■f Electronic address: |fabioC!iphys.ntu.edu.tw| 



placement fields possess logarithmic correlation functions 
at long distances, so that the memory of the crystalline 
structure is lost over large distances. In other words, the 
Bragg peaks of the world crystal are not (5-function-like, 
but display the typical behavior of a quasi-long-range or- 
der, similar to the order in a Kosterlitz-Thousless tran- 
sition in two-dimensional superfluids [23| . 

The purpose of this note is to study the generalized 
uncertainty principle (GUP) associated with the quan- 
tum physics on the world crystal and to derive physical 
consequences related to micro black hole physics. In view 
of the fact that micro black holes might be formed at en- 
ergies as low as the TeV range [1^ [l^, — which will 
be shortly available in particle accelerators such as the 
LHC, it is hoped that the presented results may be more 
than of a mere academic interest. 

The structure of our paper is as follows: In Section HIl 
we present some fundamentals of a differential calculus 
on a lattice that will be needed in the text. In Section ITTTl 
we construct position and momentum operators on a ID 
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lattice and compute their commutator. We then demon- 
strate that the usual Weyl-Heisenberg algebra Wi for p 
and X operators is on a ID lattice deformed to the Eu- 
clidean algebra E{2). By identifying the measure of un- 
certainty with a standard deviation we derive the related 
GUP on a lattice. This is done in Section IIVI There 
we focus on two critical regimes: long-wave regime and 
the regime where momenta are at the border of the first 
Brillouin zone. Interestingly enough, our GUP implies 
that quantum physics of the world-crystal universe be- 
comes "deterministic" for energies near the border of the 
Brillouin zone. In view of applications to micro black 
hole physics, we derive in Section |V] the energy-position 
GUP for a photon. Implications for micro black holes 
physics are discussed in Section IVII There we derive a 
mass-temperature relation for Schwarzschild micro black 
holes. On the phenomenological side, the latter pro- 
vides a nice resolution of a long-standing puzzle: the final 
Hawking temperature of a decaying micro black hole re- 
mains finite, in contrast to the infinite temperature of 
the standard result where Heisenberg's uncertainty prin- 
ciple operates. Besides, the final mass of the evaporation 
process is zero, thus avoiding the problems caused by the 
existence of massive black hole remnants. Entropy and 
heat capacity are discussed in Section IVIII Finally, in 
Section rVIIII we outline a connection of our results with 
't Hooft's approach to deterministic quantum mechanics 
and with deformed (or double) special relativity. Sec- 
tion |TX] is devoted to concluding remarks. For complete- 
ness, we present in Appendix an alternative derivation of 
the micro black hole mass-temperature formula. 

II. DIFFERENTIAL CALCULUS ON A 
LATTICE 

In this section we quickly review some features of a 
differential calculus on a ID lattice. An overview dis- 
cussing more asp ects of such a calculus can be found, 
e.g., in Refs. |l8l. [l9l. [20j. Independent and very elegant 
derivation of these can be also done in the framework of 
a non-commutative geometry [2§, [s^, [3l| . 

On a lattice of spacing e in one dimension, the lattice 
sites lie at x„ = ne where n runs through all integer 
numbers. There are two fundamental derivatives of a 
function f{x): 

(V/)(a:) = + , 

(V/)(x) = . (1) 

They obey the generalized Leibnitz rule 

{Wfg)(x) - {Wf){x)g(x)+f{x + e){Wg){x), 
{Wfg)(x) = {Wf){x)g{x)+f{x~e){Wg)(x). (2) 

On a lattice, integration is performed as a summation: 

/ dxfix) EE e^/(x), (3) 



where x runs over all a;„. 

For periodic functions on the lattice or for functions 
vanishing at the boundary of the world crystal, the lat- 
tice derivatives can be subjected to the lattice version of 
integration by parts: 

Y^j(x)Vg{x) = -^5(^)V/(x), (4) 

X X 

J2fi^)^9i^) = -J2gix)Vfix). (5) 

X X 

One can also define the lattice Laplacian as 
WW fix) = VV/(x) = l[f{x+e)-2f{x)+fix-e)] , (6) 

which reduces in the continuum limit to an ordinary 
Laplace operator d^. Note that the lattice Laplacian can 
also be expressed in terms of the difference of the two 
lattice derivatives: 

WW fix) = - [Wf{x) ~ Wf{x)] . (7) 
e 

The above calculus can be easily extended to any num- 
ber D of dimensions [ll,[il[2l|. 

III. POSITION AND MOMENTUM 
OPERATORS ON A LATTICE 

Consider now the quantum mechanics (QM) on a ID 
lattice in a Schrodinger-likc picture. Wave function are 
square-integrable complex functions on the lattice, where 
"integration" means here summation, and scalar prod- 
ucts are defined by 

(/Iff) = eJ2f*ix)gix). (8) 

It follows from Eq. (g]) that 

(/IV5) - -(V/I5), (9) 

so that (iV)^ — iW , and neither iW nor iW are hermitian 
operators. The lattice Laplacian ([S]), however, is hermi- 
tian. 

The position operator acting on wave functions of 
X is defined by a simple multiplication with x: 

{XJ){x) = xf{x). (10) 

Similarly we can define the lattice momentum operator 
Pc. In order to ensure hermiticity we should relate it to 
the symmetric lattice derivative [l^ [13, H^l ■ Using ^ 
we have 

{Pj){x) = A[(V/)(x) + (V/)(a;)] 

= ^Jfi^ + e)-fix-e)]. (11) 
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For small e, this reduces to the ordinary momentum op- 
erator p = —ihdx, or more precisely 



(12) 



The "canonical" commutator between and on the 
lattice reads 

([l„P,]/)(x) = |[/(:, + + 

= ih{ij){x). (13) 

The last line defines a lattice- version of the unit operator 
as the average over the two neighboring sites. Note that 
all three operators X^, Pg, and le are liermitian under the 
scalar product 

It was noted in [s^] that the operators X^, and 
ie generate the Euclidean algebra E{2) in 2D. Indeed, 
setting M = eX^, Pi = eP^/H and Pj = we obtain 

[M,Pi] = ^P2, [Af,P2] = -iPi, [Pi,P2] = 0. 

The generator M corresponds to a rotation, while Pi 
and P2 represent two translations. In the limit e — > 0, 
the Lie algebra of E{2) contracts to the standard Weyl- 

Heisenberg algebra Wi: Xe ^ x, Pe ^ p, le ^ Thus 
ordinary QM is obtained from lattice QM by a contrac- 
tion of the E{2) algebra, with the lattice spacing e playing 
the role of the deformation parameter. 

All functions on the lattice can be Fourier-decomposed 
with wave numbers in the Brillouin zone: 



dk 



27r 



/(fc)e^^- , 



with the coefficients 

/>) = e^/(x)e-^'=- 



(14) 



(15) 



This implies the good-old de Broglie relation 

mm - nfc/(fc), (16) 

and its lattice version 

(-zV/)(fc) = i^/(fc), (-zV/)(fc) = i?/(fc), (17) 
with the eigenvalues 

K = {e'^'~l)/ie = K* . (18) 
From (|17p we find the Fourier transforms of the operators 

Xe,Pe,Ie- 

(xjm = (19) 

{PJm = -^sm{ke)f{k), (20) 
{ij){k) = cos(fce)/(fc). (21) 



With the help of ([2T|) we can rewrite the commutation 
relation equivalently as 



[X„P,]f){x) = ih cos {ep/n) fix). (22) 



The latter allows to identify the lattice unit operator /g 
with cos{ep/h). Indeed, — 1 on all lattice nodes. 



IV. UNCERTAINTY RELATIONS ON LATTICE 

We are now prepared to derive the generalized uncer- 
tainty relation implied by the previous commutators. We 
shall define the uncertainty of an observable A in a state 
-tp by the standard deviation 



(AA), 



(23) 



Following the conventional Robertson-Schrodinger proce- 
dure (see, e.g., Ref. 33, 34, 35i]), we derive on the space- 
time lattice the inequality 



(AX,)^(AP,), 



> 



mx,,p,M 

{ip \ cos (ep/h) 



(24) 



For brevity we will omit in the following the subscript 
in (AA)^ and set {ip\ ■ ■ ■ IV') = (• ■ ■ 

Let us now study two critical regimes of the GUP 
the first is the long-wavelengths regime where (p)^ —> 0; 
the second regime is near the boundary of the Brillouin 
zone where (p)^ iTh/2e. To this end we first rewrite 
(cos (ep/h))^ as 



{cos (ep/h)) 



where g{p) = \ip(p)\'^. 

In the first case, g{p) is peaked around p ~ 0, so that 
the relation (j25p becomes approximately 



{cos{ep/h))^, 



2ft2 



(26) 



where p^ = {p^)ip. We should stress that expansion ((26|) 
is not an expansion in e but rather in tp/h. So if we 
speak of gip) as being peaked around p ~ we mean 
that p <C h/e. 

Applying now the identity 



we obtain from ((23 



(27) 



AXAP, > - 



1 - 



t^p'^ 



(28) 
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For mirror-symmetric states where (p)^/, = this imphes 



A.Y.AP, >|(l-|,(Ap)^ 



(29) 



Here we have substituted |...| by (...) since we assume 
that e ~ (Planckian lattice) and that Ap is close 
to zero (this is our original assumption). Therefore 
e2(Ap)V2fi" < 1. 

For Planckian lattices with the relation ^12 
neglect higher powers of e in p9|) and write 



we can 



AXeAP, > - 



(30) 



In the second case, where (p)^ — > tm /2t, i.e. near the 
border of the Brillouin zone, we use the expansion: 



(cos[7r/2 + (ep/;i-7r/2)])^ 

oo 



ri=0 



(sin(7r/2 - ep/K))^ 
dP^^P^^-^^ (2n+l)! ■ 



Under the assumption that g{p) is peaked near the border 
of the Brillouin zone, the first term in the expansion is 
dominant, and the uncertainty relation reduces to 



I - I ^P^^ 



(32) 



Since k = p/h lies always inside the Brillouin zone, we 
have (p)^ < 7rfi,/2e and can therefore in ([5^ substitute 
|...| by (...). Finally, using again p^ . we can write for 
the GUP close to the boundary of the Brillouin zone 



> i 



As the momentum reaches the boundary of the Brillouin 
zone, the right-hand sides of vanish, so that lat- 

tice quantum mechanics at short wavelengths is permit- 
ted to exhibit classical behavior — no irreducible lower 
bound for uncertainties of two complementary observ- 
ables appears! 

It is worth noting that the uncertainty relation (|33p 
leads to the same physical conclusions as those found, on 
a different ground, by Magueijo and Smolin in Ref. [36l |. 
In particular, the world-crystal universe can become "de- 
terministic" for energies near the border of the Brillouin 
zone, i.e., for Planckian energies. 

Let us remark that the scenario in which the universe 
at Planckian energies is deterministic rather than being 
dominated by tumultuous quantum fluctuations is a re- 
current theme in 't Hooft's "deterministic" quantum me- 
chanics ii, m m s m • 

It is straightforward to generalize the above formulas 
to higher dimensions. In this context, a useful inequality 
is 

ax:a|p,i > ^\mp:/\p,\) COS {eY/fi)m 

= ^meip')cos{eY/h)m, (34) 



which will be needed in the following. Here e{. . .) is the 
sign function, and 



IP. I = 



D r 



sm{e^p> /h) 



e3 



(35) 



Inequality p4p should be contrasted with inequality 
where the momentum is without an absolute value. 

In a particular case when states ip are a combination 
of only positive or only negative momentum eigenstates 
(e.g., incident or reflected particle states) we can simply 
write 

■ ;os(ey/n) m 



AXIA\PA > ^1 



1 - 2(sin2 (ey/2?i)). 



(36) 



IMPLICATIONS FOR PHOTONS 



We may now use the inequality (|36|) to derive the GUP 
for photons. 

The vector potential of a photon in the Lorentz gauge 
in 1 -|- 1 dimensions satisfies the wave equation 



1 



(37) 



A plane wave solution A^{x) = exp[i{kx — uj{k)t)] pos- 
sesses the well-known linear dispersion relation 



Lu{k) — c\k\ , 



(38) 



(33) 

with being a polarization vector. On a one- 
dimensional lattice, the operator is replaced by the 
lattice Laplacian VV, and the spectrum becomes, on ac- 
count of Eq. ([6]) and (fTT]) . 



uj{k) 



^2 [1 - cos(fce)] 



(39) 



which reduces to ([55)1 for e 0. Denoting the energy on 
the lattice hio by E^, we obtain the dispersion relation 



he 



■ (P^\ 

'in TTT 



(40) 



We can also define the associated energy operator by 
replacing p by p. 

For states V' with q{j>) sharply peaked around small p, 
we can use a spectral expansion analog of ()25p to obtain 



AE^ ~ cA|p| ~ cA|Fe| 



(41) 



Here we have neglected higher powers of momentum and 
used the fact that we deal with a Planckian lattice. In 
deriving we have also applied the cumulant expansion: 



= (2?ic/e)(|sin(pe/2?i)|),/, 



(42) 
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With the help of ([361) , (02 
long- wavelength regime 



and (|4ip we can write in the 



AXA^;, > — 



he 

2 



1 - 



2^2 c2 



{El 



(43) 



Here {E^)^ is the average quadrat of the photon energy, 
and thus the square root of it can be formally identified 
with the energy change in the detector, i.e. /^E^. From 
this follows that if the uncertainty of a photon position in 
a state ^ is AX^, then the energy of a detector changes 
at least by amount 



l\E, ~ — 



he 

2 



2h'^e^ 



{El 



1 



AX 



(44) 



per particle. Remembering the Einstein relation A£^ = 
2'Khc/X, we can interpret 47rAXg as being a lattice equiv- 
alent of photon's wavelength A. 

It is interesting to observe that in the short-wavelength 
case we can deduce from the exact GUP 



AX,A|P,| > I 



{El 



(45) 



have found a black hole on the lattice as a discretized 
version of a Schwarzschild solution. It is a pile up of 
disclinations. If the Schwarzschild radius is much larger 
than the lattice spacing e, this will not look much differ- 
ent from the well-known continuum solution. We must 
avoid too small black holes, for otherwise, completely 
new physics will set in near the center, due to the high 
concentration of defects. These will cause the "melting" 
of the world crystal at a critical defect density [s^, and 
the emerging trans-horizon general relativity would look 
completely different from Einstein's theory. 

Following the classical argument of the Heisenberg mi- 
croscope j5j| , we know that the smallest resolvable detail 
5x of an object goes roughly as the wavelength of the 
employed photons. If E is the (average) energy of the 
photons used in the microscope, then 



5x ~ — 



he 

2E 



(48) 



Conversely, with the relation one can compute the 
energy E oi a. photon with a given (average) wavelength 
A ~ 5x. As a consequence of Eq. (|43)) . we can write the 
lattice version of this standard Heisenberg formula as 



that near the border of the Brillouin zone AX^ takes the 
approximate form (cf. Eq. (201)) 



5X, 



he 
2K 



1 



2h^c^ 



(49) 



AX 



2h'^c 



{El 



In the derivation we have used that fact that 

~ 2e ■ 



AIR I < 



(46) 



(47) 



Relation represents the smallest attainable posi- 

tional uncertainty near the border of the Brillouin zone. 
It will be useful in the following two sections. 



VI. APPLICATIONS TO MICRO BLACK HOLES 

An interesting playground where one can apply the 
above lattice GUP's is the hypothetical physics of micro 
black holes. Their mass-temperature relation depends 
sensitively on the actual form of the energy-position un- 
certainty relation. From this one can deduce non-trivial 
phenomenological consequences. The passage from the 
energy-position uncertainty relation to the micro black 
hole mass-temperature relation has been intensively stud- 
ied in recent years. For definiteness we shall follow here 
the treatment of Refs. [44, 45„46„^, 48, 49, 50, 5l|. An 
alternative derivation based on the so-called Landauer 
principle will be presented in Appendix. 

We start with an assumption that the lattice spacing 
is roughly of order Planck length, i.e., e = a-^p, where 
a > is of order of unity. Let us now imagine that we 



which links the (average) wavelength of a photon to its 
energy E^. Since the lattice spacing is e = a£p and the 
Planck energy £p — hc/2tp, Eq. ([^5]) can be rewritten 
as 



5X, 



he 
2K 



a'^ipE, 
86 V 



(50) 



Let us now loosely follow the argument of Refs. [4Jj 
m, 113, m, Ei, [Ho, HH and consider an ensemble of un- 
polarized photons of Hawking radiation just outside the 
event horizon. From a geometrical point of view, it's 
easy to see that the position uncertainty of such pho- 
tons is of the order of the Schwarzschild radius Rs of the 
hole. An equivalent argument comes from considering 
the average wavelength of the Hawking radiation, which 
is of the order of the geometrical size of the hole (see 
e.g. Ref. [4§|, chapter 5). By recalling that Rs = ipiri, 
where m = M/Mp is the black hole mass in Planck units 
{Mp = Epjf?), we can estimate the photon positional 
uncertainty as 



bX^ ~ 2[lRs = 2iilpm. 



(51) 



The proportionality constant /i is of order unity and will 
be fixed shortly. According to the above arguments, m 
must be assumed to be much larger than unity, in order to 
avoid the melting transition. With (|5ip we can rephrase 
Eq. dSOl) as 



2[im, 



f 

Cyp 

e: 



8 £r> 



(52) 
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According to the equipartition principle the average en- 
ergy Eg of unpolarized photons of the Hawking radiation 
is hnked with their temperature T as 



(53) 



In order to fix fi, we go to the continuum lattice limit 
e ^ (a ^ 0), and require that formula (|52)) predicts 
the standard semiclassical Hawking temperature: 



Th = 



This fixes fi = n. 



he 



STrCksM inksRs 



(54) 



Defining the Planck temperature Tp so that £p = 
kBTp/2 and measuring all temperatures in Planck units 
'AS Q = T/Tp, we can finally cast formula ([5^ in the form 



2m = 



1 



27re 



- 2^9 



(55) 



where we have defined the deformation parameter C — 
a/(27r\/2). 

As already mentioned, in the continuum limit both 
e and a tend to zero and (|50l) reduces to the ordinary 
Heisenberg uncertainty principle. In this case Eq. ([51 
boils down to 



1 



AttS 



(56) 



This is the dimensionless version of Hawking's formula 
(151)1 for large black holes. 

Historically, the validity of ([51)1 was also postulated for 
micro black holes on the assumption that the black hole 
thermodynamics is universally valid for any black hole, 
be it formed via star collapse, or primordially via quan- 
tum fluctuations. Such an assumption is by no means 
warranted without some further input about mesoscopic 
and/or microscopic energy scales (much like in ordinary 
thermodynamics) and, in fact, we have seen that correc- 
tions should be expected at short world-crystal scales. 

It is instructive to compare our mass-temperature rela- 
tion ([55)l with the one suggested by the so-called stringy 
uncertainty relation [55!, ISa] ■ There the sign of the cor- 
rection term in (|55p is positive: 



2m = 



1 



27re 



27re . 



(57) 



The phenomenological consequences of the relation 
(155)1 are quite different from those of the stringy result 
(157)1 . In Fig. [T] we compare the two results, and add 
also the curve for the ordinary Hawking relation ()56p . 
Considering m and & as functions of time, we can fol- 
low the evolution of a micro black hole from the curves 
in Fig. 1. For the stringy GUP, the blue line predicts a 
maximum temperature 



1 

2^ 



(58a) 




0.05 0.1 



FIG. 1: Diagrams for the three mass-temperature relations, 
ours (red), Hawking's (green), and stringy GUP result (blue), 
with = y/2, as an example. As a consequence of the lattice 
uncertainty principle the evaporation ends at a finite temper- 
ature with a zero rest-mass remnant. 



and a minimum rest mass 



(58b) 



The end of the evaporation process is reached in a finite 
time, the final temperature is finite, and there is a rem- 
nant of a finite rest mass (cf. Refs. [3, HI, [H, [i^, lisl . 

From the standard Heisenberg uncertainty principle we 
find the green curve, representing the usual Hawking for- 
mula. Here the evaporation process ends, after a finite 
time, with a zero mass and a worrisome infinite tem- 
perature. In the literature, the undesired infinite final 
temperature predicted by Hawking's formula has so far 
been cured only with the help of the stringy GUP, which 
brings the final temperature to a finite value. This result 
is, however, also questionable since it implies the exis- 
tence of finite-mass remnants in the universe. Though 
by some authors such remnants are greeted as relevant 
candidates for dark matter [53], others point out that 
their existence would create further complications such 
as the entropy/information problem 58], detectability is- 
sue, or their (excessive) production in the early universe 

[illli]. 

In contrast to these results, our lattice GUP predicts 
the red curve. This yields a finite end temperature 



1 

2^ 



(59) 



with a zero-mass remnant. The mass-temperature for- 
mula ()55p thus solves at once several problems by pre- 
dicting the end of the evaporation process at a finite final 
temperature with zero-mass remnants. 

Is should be stressed that since the photon GUP (|43)) 
and (|49)l holds only for states tp where (p^)^ <C fi.^/e^, 
our reasonings are warranted only for 



Eg <C hc/e ~ £p/a 



27rC < 



1 
6 



(60) 
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This implies, in particular, that when Q is close to 0max 
our long-wavelength approximation cannot be trusted. 

To understand the behavior of the system close to 
©max we must turn to the short- wavelength limit, 
Eqs. (|33p and In this regime the momenta lie 

close to the border of the Brillouin zone {p)^ ~ {Pe)il> — 
TTh/{2e), and Eq. (gOl) iniphes 

/2 

E, ~ —he, (61) 
e 

which for the Planckian lattice, where e = a£p, gives 
~ £p/(TrQ. Considering again the uncertainty in the 
photon position as SX^^ ~ 2TT£pm (cf. Eq. (I51|)). GUP 
then predicts 

We can thus conclude that the mass of the micro black 
hole must go to zero. This is also consistent with our pre- 
vious long- wavelength considerations. The micro black 
hole therefore evaporates completely, without leaving 
remnants. 



By integrating dS we obtain S = S{Q). Just as for- 
mula ([55]) . the relation ([^ can be trusted only for 
6 ^ 6max — l/27rC. Thus, when integrating ([BS]) . we 
should do this only up to a cutoff Gmax ^ ©max- The 
additive constant in S can be then be fixed by requiring 
that S = when Gmax- This is equivalent to what 
is usually done when calculating the Hawking tempera- 
ture for a Schwarzschild black hole. There one fixes the 
additive constant in the entropy integral to be zero for 
TO = 0, so that S{m = 0) = S{Q oo) = (the mini- 
mum mass attainable in the standard Hawking effect is 
TO = 0). Thus we obtain 

where the sign was chosen in order to have a positive 
entropy. 

The integral (pS)) yields 

\ ^max / 

The entropy is always positive, and ^ for O ©max- 



VII. ENTROPY AND HEAT CAPACITY 



B. Heat Capacity 



In this section we exhibit the modified thermodynamic 
entropy and heat capacity of a black hole implied by the 
new mass-temperature formula (I55p . 



A. Entropy 

From the first law of black hole thermodynamics [60*1 
we know that the differential of the thermodynamical 
entropy of a Schwarzschild black hole reads 



dS 



dE 



(63) 



where dE is the amount of energy swallowed by a black 
hole with Hawking temperature Th- In Eq. (|63p the 
increase in the internal energy is equal to the added heat 
because a black hole makes no mechanical work when 
its entropy/surface changes (expanding surface does not 
exert any pressure). 

Rewriting Eq. (I63p with the dimensionless variables m 
and 8 we get 



dS 



kB dm 



Inserting here formula (|55p we find 



,„ _ kB_dm _ kB_ 



(64) 



de . (65) 



With entropy formulae ([65|1 and ([67|) at hand we can 
now compute the heat capacity of a (micro) black hole in 
the world-crystal. This will give us important insights on 
the final stage of the evaporation process. Again, we shall 
obtain formulae vaHd only for 6 <C ©max = l/(27r^). 

The heat capacity C of a black hole is defined via the 
relation 



dQ ^ dE = CdT . 



(68) 



The pressure exerted on the environment by the expand- 
ing black hole surface is zero. Hence we do not need to 
specify which C is meant. 

With the help of ^ and ([68]) we obtain 



C 



dT 



which yields 



C = - 



(27re)2 



(69) 



(70) 



From this clearly follows that C is always negative. 

Most condensed-matter systems have C > 0. However, 
because of instabilities induced by gravity this is gener- 
ally not the case in astrophysics [6l|, [62], especially in 
black hole physics. A Schwarzschild black hole has C < 
which indicates that the black hole becomes hotter by ra- 
diating. The result (|7D|) implies that this scenario holds 
also for micro black holes in the world-crystal. 
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In case of stringy GUP, we have to use Eq. ([57)) as the 
mass-temperature formula. The expression for the heat 
capacity then reads 



C = 



(27re)2 



(71) 



Since also here < < 8max — l/27rC, black holes have 
negative specific heat also according to the stringy GUP. 
However, stringy GUP displays a striking difference with 
respect to lattice GUP. In fact, since in principle we can 
trust Eq. ([57)1 also when Q ~ Omax = l/27rC, then from 
([7T|) we have, in such limit, C = 0. This means that 
for the stringy GUP the specific heat vanishes at the 
end point of the evaporation process in a finite time, so 
that the black hole at the end of its evolution cannot 
exchange energy with the surrounding space. In other 
words, the black hole stops to interact thermodynami- 
cally with the environment. The final stage of the Hawk- 
ing evaporation, according to the stringy GUP scenario, 
contains a Planck-size remnant with a maximal tempera- 
ture Q = Omax, but thermodynamically inert. The rem- 
nant behaves like an elementary particle — there are no 
internal degrees of freedom to excite in order to produce 
a heat absorption or emission. 

To understand the heat exchange in the last live stage 
of the world-crystal black hole we cannot use the long- 
wavelength formula (|70p . Instead, we must turn to 
Eq. Since in our scenario the micro black hole dis- 

appears at the critical temperature 0inax, it is more ap- 
propriate to write the mass-temperature formula (|46| in 
the form 



e(e,„ax - e) — c ( 1 - ^ 



(72) 



where 9{t) is the Heaviside step function. For the specific 
heat this implies 



C 



(73) 



So, in contrast to the stringy result, a world-crystal black 
hole exchanges heat with its environment by radiation 
until the last moment of its existence, and unlike the 
Schwarzschild black hole, the heat exchange with the en- 
vironment increases in the final stage of its evaporation. 
In addition, because of the 0- function in ([75]) , the transi- 
tion from the universe with the world-crystal black hole 
to the one without it is of first order. 



VIII. FURTHER APPLICATIONS 

So far we have studied the consequence of the GUP on 
the micro black holes. Let us briefly mention two further 
applications. 



A. 't Hooft's proposal 

The first application relates to 't Hooft's proposal 
which purports to justify that our quantum world is 
merely a low-energy limit of a deterministic system op- 
erating at a deeper, perhaps Planckian, level of dynam- 
ics [stI, [11]. As a deterministic substrate 't Hooft has 
proposed various cellular automata (CA) models. 

In general, a CA is an array of cells forming a discrete 
lattice. All cells are typically equivalent and can take one 
of a finite number of possible discrete states. Like space, 
time is discrete as well. At each time step every cell up- 
dates its state according to a transition rule which takes 
into account the previous states of cells in the neighbor- 
hood, including its own state. In this sense the evolution 
is deterministic. 

One of the simplest CA considered by 't Hooft is the 1- 
dimensional periodic CA with 4-state cells, and with the 
nearest neighbor (N-N) transition rule, see Fig. O This 







-^3 



G &- 







&- 



-O 



FIG. 2: CA with discrete time evolution described by 
Eq. (|74|l and with the periodicity condition ai = cri+4. The 
right-hand shows an equivalent graphical representation. 



"clock like" CA can be generally described with 2N + 1 
cells (Tj {i ~ —N,..., N) each with two possible states 
{0,1} (cell is, e.g., white or black). The discrete time 
evolution with the elementary time step St is described 
by the N-N transition rule 

(o-i-i, o-j, (jj+i) ^ o-^ = cri{t + St): 

(0,0,0)^0, (0,0,l)-.0, (0,l,0)-.0, 
(1,0,0)^1, (1,1,0)^0, (0,1,1)^0, 
(1,0,1)^0, (1,1,1)^0, (74) 

with the Born-von Karman periodicity condition (Ji = 

The cells can be algebraically represented by orthonor- 
mal vectors 







(7-N 







UN 



viy vo/ voy 

d-N — cttv+i- On the basis spanned by <Ji the 
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elementary-time step evolution operator is: 



U{dt = r) = e 



-iHSt 



— g 2N + 1 



1 



V 



(75) 



1 Oj 



which, among others, defines the Hamiltonian H. The 
pre-factor e~'2jv+i jg known as 't Hooft's phase conven- 
tion. Because [/2Af+i _ ^ ^^^g ^g^^^ diagonalize [/ as 



e *2jv+i diagle '^w+i ^ 



1, 



If we denote the eigenstates of -ff as we find that 

AT 

ae , (76) 



V2N 



1 ^ 

exp 



i=-N 



2N+1 



with n — —N, . . . , N and the ensuing "energy" spectrum 
reads 

27r 



(2iV4- l)(5i 



(77) 



The energy values i?„ resemble the spectrum of the har- 
monic oscillator, except that the n's are bounded and can 
attain negative values. We shall be coming back to this 
issue shortly. 

Position and momentum can be represented by opera- 
tors with the matrix representations 

/{-N + l)e ••• \ 

{-N + 2)e 



••• A^e 

••• -Ne/ 

/ 1 

-10 1 



n = — 



1 

2ie 



0-101 

• • • 
V 1 ••• 



-1\ 






1 
-1 / 



(78) 



where e — 27r/ (2iV -I- 1). With these we obtain the com- 
mutator 

/O 1 ••• 1\ 
10 1 ••• 
10 1 ••• 



••• 

Vi ... 



1 

1 0/ 



(79) 



In deriving (|79p we have used the the periodicity condi- 
tion a-i — <72N- By comparing this with the evolution 
operator ([75l) we have 



(80) 



[X„P,] = icos[{H -L0N/2)St], 



which for small St gives 



1 



5f 



n) 



In addition, we deduce from (|75p and ([78|l that 

Bin[{H ~ujN/2)5t] = eP, , (82) 



which is compatible with the result l|20|) . From ((82)) fol- 
lows that H depends only on but not on X^. So 
and H are simultaneously diagonalizable. By defining 
the operators K± as 



K+ = 
(k^ = X|), so that 
2iV + 1 



P. 



Pee 



(83) 



K+\n) 
k\n) 



2tt 
2N + 1 



2tt 



2tt 



2N+1 

2tt 
2N+1 



(n-l) 



(84) 



one can persuade itself that H and K± close the deformed 
algebra which in the large N limit (i.e., in the small e or 
St limit) reduces to 



[K+,K^] = -— , (85) 

UJ 



with uj = Woo. 

Note that for large N one can identify ([77]) . (l84|) . and 
([85|l with the representation of SU(1,1) known as the 
discrete series D^^® D~^., (cf. Ref. [6^). Generally, the 

Lie algebra ® is defined through the relations: 

L3\k,m) = {m + k)\k,m) , 
L+\k,m) = v^(m + 2k){m + 1) |fc, m -|- 1) , 
L_|fc,TO) = V(n + 2fc 
[L3,i±] = ±^±, 



1)to \k, m — 1) , 

-2i3 . (86) 



Here m + k = ±fc, ±(fc 1), ±{k + 2),... and k = 
^ 1 - is the so-called Bargmann index which labels 



1 - 2 

2 ' ' 2 ' ' ■ 



1/2 



(87) 



the representations. From this we have that D^^^ ( 
corresponds to 

Lslhm) - (m + l/2)|i,m), 
i+|i,m) = (m + l)|i,m + l), 
L_|i,m) = m|i,m- 1) , 



Identification with the large- limit of ([77]) and is 
established when we identify in ([87)l with 07X3, if± 
with L-t, and set m = n. 

At this stage one can invoke 't Hooft's loss of informa- 
tion condition [13, HI], and project out the negative part 
of the spectra. A plausible rationale for this step can 
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be found, e.g., in irreversibility of computational process 
due to a finite storage capacity |6J] . 

After the negative energy spectrum is removed 
(erased), we obtain only the positive discrete series 
£'^^2 representation where m = 0,1,2,...), and 

the Hamiltonian morphs into a non-negative spectrum 
Hamiltonian H~^. 

The usual H^(l)-algebra of the quantum harmonic os- 
cillator emerges after we introduce the following mapping 
in the universal enveloping algebra of SU{1, 1): 



1 



1/2 



it = L+ 



1 



1/2 



(88) 



The latter gives a one-to-one (non-linear) mapping be- 
tween the deterministic cellular automaton system (with 
information loss) and the quantum harmonic oscillator. 
The reader will recognize the mapping Eq. (j88p as 
the non-compact analog 65] of the well-known Holstein- 
Primakoff representation for SU{2) spin systems [66i] . 

Our operators X^, and may be viewed as the 
same cellular automaton E{2) algebra as discussed in 
Section Hill Thus we can conclude that in the Planckian 
scale the system must behave deterministically — which 
is one of the defining property of cellular automata. It is 
only at low energies when the loss of information leads 
to the emergent degrees of freedom resulting in the usual 
quantum mechanical description 



B. Double special relativity 

The second application relates to the idea of double (or 
doubly or deformed) special relativity (DSR) (see, e.g., 
Refs. m, [G^I). The general idea is that if the Planck 
length is a truly universal quantity, then it should look 
the same to any inertial observer. This demands a mod- 
ification (deformation) of the Lorenz transformations, to 
accommodate an invariant length scale. In Ref. [36| the 
nonlinearity of the deformed Lorenz transformations lead 
the authors to novel commutators between spacetime co- 
ordinates and momenta, depending on the energy 



E 

y 



(89) 



where E is the energy scale of the particle to which the 
deformed Lorenz boost is to be applied, while £p is the 
Planck energy. This suggests that they have an energy- 
dependent Planck "constant" h{E) ~ h{\ — E/Ep). Their 
model also implies that h{E) — + for ^ £p. For en- 
ergies much below that Planck regime, the usual Hcisen- 
berg commutators are recovered, but when E ~ £p one 
has h{£p) ~ 0. So the Planck energy is not only an in- 
variant in this model, but the world looks also apparently 
classical at the Planck scale, similarly as in 't Hooft's pro- 
posal. 



The connections of the DSR model with our proposal 
are at this point self evident. Our GUP (HH), ([Ml) im- 
plies that, at the boundary of the Brillouin zone, when 
{p)tl) ~^ h7T/2e, i.e. for Planck energies 
the fundamental commutator vanishes 



and since 



AXAP, > 0, 



(90) 



(91) 



lattice quantum mechanics at short wavelengths allows 
for classical behavior, that is uncertainties of two com- 
plementary observables can be simultaneously zero. 

However, if we express the fundamental commutator 
P2|) of our model in terms of energy, using the exact 
relation (|40|) . we find (for e — a£p) 



fix) ^ ihli 



(92) 



This means that the deforming term in our model is 
quadratic in the energy, instead of the linear dependence 
in the energy of the DSR model ([89|. 



IX. DISCUSSION AND SUMMARY 

It should be noted that the present lattice generaliza- 
tion of the uncertainty principle is not an approximate 
description, but it is an exact formula necessarily implied 
by our model of lattice space time. The great majority 
of the GUP research has always borrowed the deformed 
commutator [x,p] = ih{l + Kp^) either from strin g th eor y , 
or from heuristic arguments about black holes [53, \EM ■ 
To be precise, even in string theory the formula ex- 
pressing the GUP is not derived from the basic features 
of the model, but instead it is deduced from high-energy 
gedanken experiments of string scatterings. In contrast 
to this we have derived all results from a simple lattice 
model of spacetime, and from the analytic structure of 
the basic commutator ([22]). 

We have calculated the uncertainties on a crystal-like 
universe whose lattice spacing is of the order of Planck 
length — the so-called world crystal. When the ener- 
gies lie near the border of the Brillouin zone, i.e., for 
Planckian energies, the uncertainty relations for position 
and momenta do not pose any lower bound on the as- 
sociated uncertainties. Hence the world crystal universe 
can become "deterministic" at Planckian energies. In this 
high-energy regime, our lattice uncertainty relations re- 
semble the double special relativity result of Magueijo 
and Smolin. 

The scenario in which the universe at Planckian en- 
ergies is deterministic rather than being dominated by 
quantum fluctuations is a starting point in 't Hooft's "de- 
terministic" quantum mechanics. 
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With the generalized uncertainty relation at hand we 
have been able to derive a new mass-temperature relation 
for Schwarzschild micro black holes. In contrast to stan- 
dard results based on Heisenberg or stringy uncertainty 
relations, our mass-temperature formula predicts both 
finite Hawking's temperature and a zero rest-mass rem- 
nant at the end of the evaporation process. Especially the 
absence of remnants is a welcome bonus which allows to 
avoid such conceptual difficulties as entropy/information 
problem or why we do not experimentally observe the 
remnants that must have been prodigiously produced in 
the early universe. 
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Apart from the mass-temperature relation we have also 
computed two relevant thermodynamic characteristics, 
namely entropy and heat capacity. Particularly the heat 
capacity provided an important insight into the last life 
stage of the world-crystal micro black holes. In contrast 
to the stringy result, our result indicates that world- 
crystal micro black hole exchanges heat with its environ- 
ment (radiate) till the last moment of its existence, and 
unlike the Schwarzschild micro black hole, the heat ex- 
change with the environment increases in the final stage 
of the evaporation. In addition, the transition from the 
universe with the world-crystal micro black hole to the 
one without is of the first order. 



Since the world crystal physics allows for deterministic 
description of the physics at Planckian energies, we have 
included in this paper a discussion of 't Hooft's periodic 
automaton model which gives at low energy scales rise 
to a genuine quantum harmonic oscillator. In addition, 
such an automaton has a close connection with our world 
crystal paradigm. Here we have re-derived 't Hooft's re- 
sult in a new way. In contrast to 't Hooft derivation [ssj 
we have matched the algebra of the automaton variables 
with the SU {1,1) algebra, and in contrast to Ref. [6^ 
we have worked with a different set of dynamical vari- 
ables. In the contraction limit (i.e., in the limit of many 
lattice spacings — low energy limit) we have recovered 
the canonical W^(l)-algebra and were able to identify the 
"emergent" harmonic oscillator variables. 



There are several aspects of the double special relativ- 
ity that are worth noting in the connection with our gen- 
eralized uncertainty relation. Essentially, we have seen 
that the fundamental commutator in DSR as well as in 
our lattice GUP goes to zero at Planck energy. For both 
models, the world should therefore be manifestly "classi- 
cal" in the Planck regime, a feature very different from 
the common believe. This is a striking prediction sup- 
ported by both models, although the lines of thought 
followed in the two research paths are completely differ- 
ent and independent. Moreover, this aspect presents also 
a strong resemblance with the results obtained in the re- 
search line of "deterministic" quantum mechanics. 



Appendix: Landauer principle 

Here we wish to provide an alternative derivation of 
the mass-temperature formula (j55p based on Landauer 's 
principle. To this end we consider an ensemble of un- 
polarized photons that are going to deliver to a micro 
black hole one single bit of information per particle. In 
order to be sure that each photon delivers only one bit of 
information — namely the information that it is there, 
somewhere inside in the black hole, its position uncer- 
tainty must be "maximal", i.e., it should not be smaller 
than Schwarzschild's radius Rs as otherwise the photon 
would deliver to black hole also extra bits of information 
concerning its entry point (or better sector) on the hori- 
zon. At the same time its wavelength should not be big- 
ger than Rs, as otherwise the photon would bounced off 
the black hole without getting trapped. In this view, the 
position uncertainty of a photon in the ensemble must be 
of order of Schwarzschild's radius Rs i.e., AX^ ~ 
Factor /x is Bekenstein's deficit coefficient which ensures 
a correct Hawking's formula in a continuum limit. 

An extra bit of information added to the micro black 
hole will increase its energy at least by amount AE^ so 
that (cf. dS])) 



he 

2 



1 



{AE^f 



(93) 



In the following we denote AE^ simply as E^ to stress 
that AEc an energy increas due to one photon. With the 
explicit form for Planck's energy 



£:„ = 4r « O.ei-lOi^GeV, 



the relation can be cast to 

he a'^e„E, 



AX, 



2Ef 8£rt 



(94) 



(95) 



If we use further the fact that, Rs = ip m, where m is 
the relative mass of the black hole in Planck units, i.e., 
m — M/Mp {Mp — £p/c^), we can rewrite (1951) as 



2m/i 



£jy a E^ 



Ef S>£r, 



(96) 
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According to the Landauer principle [g^ , when a single 
bit of information is erased (like in the black hole) the 
amount of energy dissipated into environment is at least 
kBT\ii2, where ks is Boltzmann's constant and T is the 
temperature of the erasing environment (in our case the 
micro black hole). Since the liberated energy per bit of 
lost information can not be grater the energy of the 
carrier photon we have that 

E, ~ keT. (97) 

Relation ([W]) basically expresses equipartition law for 



an unpolarized photon in the outgoing Hawking radia- 
tion. Defining the Planck temperature Tp = 2£p/kB « 
3 • lO'^^ K, and measuring the temperature in terms of 
Planck units as a relative temperature — T/Tp, we 
can rewrite Eq. (j96[) in the form 

2m = ^-2<2e. (98) 

where we identify ^ — a/(2V27r) and set /i = tt, in order 
to agree with (|55p and with Hawking's formula (j56[) in 
the continuum limit. 
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